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1. Introduction 

Numerical computations are still very important in computer applica- 
tions. But until recently there has been a discrepancy between numer- 
ical methods and software/hardware tools for scientific calculations. 
In particular, numerical programming was not much influenced by 
the progress in Mathematics, programming languages and technology. 
Modern tools for numerical calculations are not unified, standardized 
and reliable enough. It is difficult to ensure the necessary accuracy and 
safety of calculations without loss of the efficiency and speed of data 
processing. It is difficult to get correct and exact estimations of calcu- 
lation errors. For example, standard methods of interval arithmetic M 



do not allow to take into account error auto-correction effects IS] and, 
as a result, to estimate calculation errors accurately. 

However, new ideas in Mathematics and Computer Science lead to a 
very promising approach (initially presented in |20[|-||22[|). An essential 
aspect of this approach is developing a system of algorithms, utilities 
and programs based on a new mathematical calculus which is called 
Idempotent Analysis, Idempotent Calculus, or Idempotent Mathemat- 
ics. For many problems in optimization and mathematical modeling 

* The work was supported by the RFBR grant 99-01-01198. 
Submitted to Reliable Computing. 
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this Idempotent Analysis plays the same unifying role as Functional 
Analysis in Mathematical Physics, see, e.g., ]14J], [O], p9[-pl]1 and 
surveys [15], [21|. 

Idempotent Analysis is based on replacing the usual arithmetic oper- 
ations by a new set of basic operations (such as maximum or minimum) . 
There is a lot of such new arithmetics, which are associated with 
sufficiently rich algebraic structures called idempotent semirings. It is 
very important that many problems that are nonlinear in the usual 
sense become linear with respect to an appropriate new arithmetic, i.e., 
linear over a suitable semiring (the so-called idempotent superposition 
principle 27], [28], which is a natural analog of the well-known 
superposition principle in Quantum Mechanics). This 'linearity' con- 
siderably simplifies explicit construction of their solutions. Examples 
are the Bellman equation and its generalizations, the Hamilton— Jacobi 
equation, etc. The idempotent analysis serves as a powerful heuris- 
tic tool to construct new algorithms and apply unexpected analogies 
and ideas borrowed, e.g., from mathematical physics and quantum 
mechanics. 

The abstract theory is well advanced and includes, in particular, a 
new integration theory, linear algebra and spectral theory, idempotent 
functional analysis, idempotent Fourier transforms, etc. Its applications 
include various optimization problems such as multi-criteria decision 
making, optimization on graphs, discrete optimization with a large 
parameter (asymptotic problems), optimal design of computer systems 
and computer media, optimal organization of parallel data processing, 
dynamic programming, applications to differential equations, numerical 
analysis, discrete event systems, computer science, discrete mathemat- 
ics, mathematical logic, etc. (see, e.g., [|, [§, [§-(II|, 0, gj], [§7||-|11 
and references therein). 

It is possible to implement this new approach to scientific and nu- 
meric calculations in the form of a powerful software system based on 
unified algorithms. This approach ensures the arbitrary necessary ac- 
curacy and safety of numerical calculations and working time reduction 
for developers and testers of algorithms because of software unification. 

Our approach uses techniques of object oriented and functional pro- 
gramming (see, e.g., |16[| ), which are very convenient for the design 
of our (suggested) software system. Computer algebra techniques j| 
are also used. Modern techniques of systolic processors and VLSI real- 
izations of numerical algorithms including parallel algorithms of linear 
algebra (see, e.g., [l^], |32]]) are convenient for effective implementations 
of the proposed approach to hardware design. 

There is a regular method based on the idempotent theory for con- 
structing back-end processors and technical devices intended for real- 
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ization of basic algorithms of idempotent calculus and mathematics 
of semirings. These hardware facilities can increase the speed of data 
processing. 

2. Mathematical objects and their computer representations 

Numerical algorithms are combinations of basic operations. Usually 
these basic operations deal with 'numbers'. In fact these 'numbers' 
are thought of as members of some numerical domains (real numbers, 
integers etc.). But every computer calculation deals with finite models 
(or finite computer representations) of these numerical domains. For 
example, integers can be modeled by integers modulo 2 n , real numbers 
can be represented by rational numbers or floating-point numbers etc. 
Discrepancies between mathematical objects (e.g., 'ideal' numbers) and 
their computer models (representations) lead to calculation errors. 

Due to imprecision of sources of input data in real-world problems, 
the data usually come in the form of confidence intervals or other num- 
ber sets rather than exact quantities. Interval Analysis (see, e.g., 0|) 
extends operations of traditional calculus from numbers to number in- 
tervals to make possible processing such imprecise data and controlling 
rounding errors in computational mathematics. 

However, there is no universal model that would be good in all 
cases, so we have to use varieties of computer models. For example, 
real numbers can be represented by the following computer models: 

— standard floating-point numbers, 

— double precision floating-point numbers, 

— arbitrary precision floating-point numbers, 

— rational numbers, 

— finite precision rational numbers, 

— floating-slash and fixed-slash rational numbers, 

— interval numbers, 

— and others. 

When examining an algorithm it is often useful to have a possibility 
to change computer representations of input/output data. For this aim 
the corresponding algorithm (and its software implementation) must 
be universal enough. 
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3. Universal algorithms 

It is very important that many algorithms do not depend on particular 
models of a numerical domain and even on the domain itself. Algo- 
rithms of linear algebra (matrix multiplication, Gauss elimination etc.) 
are good examples of algorithms of this type. 

Of course, one algorithm may be more universal than another al- 
gorithm solving the same problem. For example, numerical integration 
algorithms based on the Gauss-Jacobi quadrature formulas actually 
depend on computer models because they use finite precision constants. 
On the contrary, the rectangular formula and the trapezoid rule do not 
depend on models and in principle can be used even in the case of 
idempotent integration (see below). 

The so-called object oriented software tools and programming lan- 
guages (like C ++ and Java, see, e.g., |p6| ) are very convenient for 
computer implementation of universal algorithms. 

In fact there are no reasons to restrict ourselves to numerical do- 
mains only. Actually it may be a ring of polynomials, a field of rational 
functions, or an idempotent semiring. The case of idempotent semirings 
is extremely important because of numerous applications. 

4. Idempotent correspondence principle 

There is a nontrivial analogy between Mathematics of semirings and 
Quantum Mechanics. For example, the field of real numbers can be 
treated as a 'quantum object' with respect to idempotent semirings, 
which in turn can be treated as 'classical' or 'semi-classical' objects 
with respect to the former. 

Let R be the field of real numbers and R+ the subset of all non- 
negative numbers. Consider the following change of variables: 

u I— ► w = h In u, 

where u € R+ \ {0}, h > 0; thus u = e w / h , w G R. Denote by the 
additional element — oo and by S the extended real line R U {0}. The 
above change of variables has a natural extension to the whole S 
by D h (Q) = 0; also, we denote D h (l) = = 1. 

Denote by Sh the set S equipped with the two operations ©j, (gener- 
alized addition) and 0^ (generalized multiplication) such that is a 
homomorphism of {R+, +, •} to {S, ffi^, Qh}- This means that Dh{u\ + 
U2) = D h {u\) ® h D h (u 2 ) and D h (u\ ■ u 2 ) = D h {u\) Q h D h (u 2 ), i.e., 
wi Qh w 2 = w\ + W2 and W\ ©^ w 2 = h\n(e wi ^ h + e W2 ^ h ). It is easy to 
prove that w\ ©^ w 2 — > maxjwi, ^2} as h — > 0. 
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Denote by R max the set S = R U {0} equipped with operations 
© = max and = +, where = — oo, 1 = as above. Algebraic 
structures in R + and Sh are isomorphic; therefore R max is a result of a 
deformation of the structure in R + . 

We stress the obvious analogy with the quantization procedure, 
where h is the analog of the Planck constant. In these terms, R+ (or 
R) plays the part of a 'quantum object' while R m ax acts as a 'classical' 
or 'semi-classical' object that arises as the result of a dequantization of 
this quantum object. 

Likewise, denote by R m i n the set R U {0} equipped with operations 
© = min and = +, where = +oo and 1 = 0. Clearly, the 
corresponding dequantization procedure is generated by the change of 
variables u ^ w = — hhiu. 

Consider also the set R U {0, 1}, where = — oo, 1 = +oo, together 
with the operations © = max and = min. Obviously, it can be 
obtained as a result of a 'second dequantization' of R or R + . 

The algebras presented in this section are the most important exam- 
ples of idempotent semirings, the central algebraic structure of Idem- 
potent Analysis. 

Consider a set S equipped with two algebraic operations: addition 
© and multiplication 0. The triple {S, ©, 0} is an idempotent semiring 
if it satisfies the following conditions (here and below, the symbol * 
denotes any of the two operations ©, 0): 

— the addition © and the multiplication are associative: x*(y-kz) = 
(x *y) * z for all x,y, z G S; 

— the addition © is commutative: x © y = y © x for all x, y G S; 

— the addition © is idempotent: x © x = x for all x G S; 

— the multiplication is distributive with respect to the addition ©: 

x (y © z) = (x y) © (x z) and (x © y) z = (x z) © (y z) 
for all x,y,z G S. 

A unity of a semiring S is an element 1 G S such that for all x G S 

101=101=1. 

A zero of a semiring S is an element G S such that 0/1 and for 
all x G S 

0®x = x, O0x = x0O = O. 

A semiring S is said to be commutative if x y = y x for all 
x, y G S. A commutative semiring is called a semifield if every nonzero 
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element of this semiring is invertible. It is clear that K max and M m i n are 
semifields. 

Note that different versions of this axiomatics are used, see, e.g., 

i, a, @, §, m, 0-i 



17 1, [21], [31 1 and some literature 
indicated in these books and papers. Many nontrivial examples of idem- 
potent semirings can be found, e.g., in ||, Q, Jl2 |, Q, |14], p7 ], 
JH|], [23], p5| ], [31]. For example, every vector lattice or ordered group 
can be treated as an idempotent semifield. 

The addition © defines the following standard partial order on S: 
x :f> y if and only if x © y = y. If S contains a zero 0, then ^ x for all 
x G S. The operations © and are consistent with this order in the 
following sense: if x ^ y, then x * z ^ y * z and z * x H z * y for all 

The basic object of the traditional calculus is a function defined on 
some set X and taking its values in the field M (or C); its idempotent 
analog is a map X — > S, where X is some set and S = M m i n , M m a X! 
or another idempotent semiring. Let us show that redefinition of basic 
constructions of traditional calculus in terms of Idempotent Mathemat- 



ics can yield interesting and nontrivial results (see, e.g., [17], [21], [23], 
[p5|[, for details, additional examples and generalizations). 

Example 1. Integration and measures. To define an idempo- 
tent analog of the Riemann integral, consider a Riemann sum for a 
function ip(x), x E X = [a, b], and substitute semiring operations © 
and for operations + and • (usual addition and multiplication) in its 
expression (for the sake of being definite, consider the semiring R max ): 



N 



N 



y~] (pjxj) ■ Aj i-> ffiy?^) Aj = max (v?(xi) + Aj), 

z — • i=\,....N 
i=l i=l 

where a = x$ < x\ < • ■ ■ < xn = b, A, 
maxj Aj — > 0, the integral sum tends to 



I.. 



N. As 



for any function <p: X 
X the set function 



(f(x) dx = sup <p(x) 

A ' x&X 

nax that is bounded. In general, for any set 



m^(B) = sup</?(x), 8cl, 

is called an M max - measure on A"; since m !fi (\J a B a ) = sup^ m v (B a ), this 
measure is completely additive. An idempotent integral with respect to 
this measure is defined as 



ip(x) dm^ 



'x 



ip(x) ip(x) dx = sup (ip(x) + <p(x)). 
x xex 
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Using the standard partial order it is possible to generalize these 
definitions for the case of arbitrary idempotent semirings. 

Example 2. Fourier-Legendre transform. Consider the topo- 
logical group G = W 1 . The usual Fourier-Laplace transform is defined 
as 

ip(x) i — ► = / e^' x ip(x) dx, 
Jg 

where exp(i£ • x) is a character of the group G, i.e., a solution of the 
following functional equation: 

f(x + y)=f(x)f(y). 

The idempotent analog of this equation is 

f(x + y) = f(x)Qf(y) = f(x) + f(y). 

Hence 'idempotent characters' of the group G are linear functions of 
the form x i— > £ ■ x = £1X1 + • • • + £, n %n- Thus the Fourier-Laplace 
transform turns into 

<p(x) i — > = ^ ■ x ip(x) dx = sup (£ • x + 

•/(.• x&G 

This is the well-known Legendre (or Fenchel) transform. 

These examples suggest the following formulation of the idempotent 



correspondence principle |f2Q] , [21]: 



There exists a heuristic correspondence between interesting, useful 
and important constructions and results over the field of real (or 
complex) numbers and similar constructions and results over idem- 
potent semirings in the spirit of N. Bohr's correspondence principle 
in Quantum Mechanics. 

So Idempotent Mathematics can be treated as a 'classical shadow 
(or counterpart)' of the traditional Mathematics over fields. 

In particular, an idempotent version of Interval Analysis can be 
constructed p5[| . The idempotent interval arithmetic appear to be re- 
markably simpler than its traditional analog. For example, in the tradi- 
tional interval arithmetic multiplication of intervals is not distributive 
with respect to interval addition, while idempotent interval arithmetics 
conserve distributivity. Idempotent interval arithmetics are useful for 
reliable computing. 



5. Idempotent linearity 

Let S be a commutative idempotent semiring. 
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The following example of a noncommutative idempotent semiring is 
very important. 

Example 3. Let Mat n (S) be a set of all S-valued matrices, i.e., 
coefficients of these matrices are elements of S. Define the sum © of 
matrices A = \\a,ij\\, B = G Mat n (S) as A © B = \\aij © bij\\ G 

Mat n (S). The product of two matrices A G Mat ra (S) and B G Mat n (S) 
is a matrix AB G M&t n (S) such that AB = | ©^Li«ife The 
set Mat n (S) of square matrices is an idempotent semiring with respect 
to these operations. If is the zero of S, then the matrix O = \\oij\\, 
where Ojj = 0, is the zero of Mat n (S); if 1 is the unity of S, then the 
matrix E = \\5ij\\, where 5ij = 1 if i = j and Sij = otherwise, is the 
unity of Mat n (S). 

Now we discuss an idempotent analog of a linear space. A set V is 
called a semimodule over S (or an S'-semimodule) if it is equipped with 
an idempotent commutative associative addition operation ®v an d a 
multiplication Qy- S x V — > V satisfying the following conditions: for 
all A, fi € S, v, w € V 

- (A© n) Q V v = \Q V (fiQ v v); 

- A Qy (v ©y w) = (A Qy v) ©y (A ©y w); 

- (A ffi/u) Q V V = (XQ V v) ©y QuQy v). 

An S'-semimodule V is called a semimodule with zero if € S and there 
exists a zero element 0y € V such that for all v € V, A € 5 

- 0y ©y f = u; 

- AQyOy = 0Q V V = Oy. 

Example 4. Finitely generated free semimodule. The sim- 
plest S'-semimodule is the direct product S n = {(ai,...,a n ) | aj € 
S, j = l,...,n}. The set of all endomorphisms S n — > S n coincides 
with the semiring Mat n (S) of all S- valued matrices (see Example 3). 

The theory of S-valued matrices, similar to the well-known Perron- 
Frobenius theory of nonnegative matrices, is well advanced and has 
very many applications, see, e.g., @, §|, 0, H, ||, Q, 



-@)- 

Example 5. Function spaces. An idempotent function space ^(X; S) 
consists of functionals defined on a set X and taking their values in an 
idempotent semiring S. It is a subset of the set of all maps X — > S such 
that if /(x), g(x) G J^(A; S) and c £ S, then (f®g)(x) = f(x)@g{x) G 
T{X;S) and (c /)(x) = c /(x) G !F(X;S); in other words, an 
idempotent function space is another example of an S-semimodule. If 
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the semiring S contains a zero element and J~{X; S) contains the zero 
constant function o(x) = 0, then the function space T{X; S) has the 
structure of a semimodule with the zero o(x) over the semiring S. If 
the set X is finite we get the previous example. 

Recall that the idempotent addition defines a standard partial order 
in S. An important example of an idempotent functional space is the 
space B{X; S) of all functions X — > S bounded from above with respect 
to the partial order ^ in S. There are many interesting spaces of this 
type including C(X;S) (a space of continuous functions defined on a 



topological space X), analogs of the Sobolev spaces, etc. (see, e.g., |17] 



[glj], |2J, ||9J-|§1| for details). 

According to the correspondence principle, many important con- 
cepts, ideas and results can be converted from usual Functional Analy- 
sis to Idempotent Analysis. For example, an idempotent scalar product 
in B(X; S) can be defined by the formula 

r® 

(tp,ip) = / tp(x) © ip(x) dx, 
Jx 

where the integral is defined as the 'sup' operation (see example 1). 

Example 6. Integral operators. It is natural to construct idem- 
potent analogs of integral operators of the form 

K : ip(y) i ^ (K(p)(x) = K(x, y) ip(y) dy, 

where tp(y) is an element of a functional space ^(Y; S), [Kip)(x) be- 
longs to a space Ti{X\ S) and K(x, y) is a function X x Y —> S. Such 
operators are linear, i.e., they are homomorphisms of the corresponding 
functional semimodules. If S = M max , then this definition turns into the 
formula 

(Ktp)(x) = sup (K(x, y) + tp(y)). 
Formulas of this type are standard in optimization theory. 



6. Superposition principle 

In Quantum Mechanics the superposition principle means that the 
Schrodinger equation (which is basic for the theory) is linear. Simi- 
larly in Idempotent Mathematics the idempotent superposition prin- 
ciple means that some important and basic problems and equations 
(e.g., optimization problems, the Bellman equation and its versions 
and generalizations, the Hamilton- Jacobi equation) that are nonlinear 
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in the usual sense can be treated as linear over appropriate idempotent 
semirings, see 0- 

Linearity of the Hamilton-Jacobi equation over M m i n (and K max ) 
can be deduced from the usual linearity (over C) of the correspond- 
ing Schrddinger equation by means of the dequantization procedure 
described above (in Section 4). In this case the parameter h of this 
dequantization coincides with ih , where h is the Planck constant; so 
in this case h must take imaginary values (because h > 0; see |23| for 
details). Of course, this is closely related to variational principles of 
mechanics. 

The situation is similar for the differential Bellman equation, see 
0- 

It is well-known that discrete versions of the Bellman equation can 
be treated as linear over appropriate idempotent semirings. The so- 
called generalized stationary (finite dimensional) Bellman equation has 
the form 

X = AX © B, 

where X, A, B are matrices with elements from an idempotent semiring 
and the corresponding matrix operations are described in example 3 
above; the matrices A and B are given (specified) and X is unknown. 

B.A. Carre H used the idempotent linear algebra to show that differ- 
ent optimization problems for finite graphs can be formulated in unified 
manner and reduced to solving these Bellman equations, i.e., systems 
of linear algebraic equations over idempotent semirings. For example, 
Bellman's method of solving shortest path problems corresponds to 
a version of Jacobi's method for solving systems of linear equations, 
whereas Ford's algorithm corresponds to a version of Gauss-Seidel's 
method. 



7. Correspondence principle for computations 



Of course, the idempotent correspondence principle is valid for algo- 
rithms as well as for their software and hardware implementations 



221. Thus: 



// we have an important and interesting numerical algorithm, then 
there is a good chance that its semiring analogs are important and 
interesting as well. 

In particular, according to the superposition principle, analogs of 
linear algebra algorithms are especially important. Note that numer- 
ical algorithms for standard infinite-dimensional linear problems over 
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idempotent semirings (i.e., for problems related to idempotent inte- 
gration, integral operators and transformations, the Hamilton-Jacobi 
and generalized Bellman equations) deal with the corresponding finite- 
dimensional (or finite) 'linear approximations'. Nonlinear algorithms 
often can be approximated by linear ones. Thus the idempotent linear 
algebra is a basis for the idempotent numerical analysis. 

Moreover, it is well-known that linear algebra algorithms are conve- 
nient for parallel computations; their idempotent analogs admit paral- 
lelization as well. Thus we obtain a systematic way of applying parallel 
computation to optimization problems. 

Basic algorithms of linear algebra (such as inner product of two vec- 
tors, matrix addition and multiplication, etc.) often do not depend on 
concrete semirings, as well as on the nature of domains containing the 
elements of vectors and matrices. Algorithms to construct the closure 
A* = 1 ^4 A 2 © ■ ■ ■ © A n © • • • = @™ =l A n of an idempotent matrix 
A can be derived from standard methods for calculating (1 — A)^ 1 . 
For the Gauss-Jordan elimination method (via LU-decomposition) this 
trick was used in []32| , and the corresponding algorithm is universal and 
can be applied both to the Bellman equation and to computing the 
inverse of a real (or complex) matrix (1 — A). Computation of A" 1 can 
be derived from this universal algorithm with some obvious cosmetic 
transformations . 

Thus it seems reasonable to develop universal algorithms that can 
deal equally well with initial data of different domains sharing the same 
basic structure pi] ], p2[ . 



8. Correspondence principle for hardware design 

A systematic application of the correspondence principle to computer 
calculations leads to a unifying approach to software and hardware 
design. 

The most important and standard numerical algorithms have many 
hardware realizations in the form of technical devices or special proces- 
sors. These devices often can be used as prototypes for new hardware 
units generated by substitution of the usual arithmetic operations for its 
semiring analogs and by adding a representation of neutral elements 
and 1 (the latter usually is not difficult). Of course, the case of 
numerical semirings consisting of real numbers (maybe except neu- 
tral elements) is the most simple and natural p0[|-[p2|. Note that for 
semifields (including R max an d M mm ) the operation of division is also 
defined. 
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Good and efficient technical ideas and decisions can be transposed 
from prototypes into new hardware units. Thus the correspondence 
principle generates a regular heuristic method for hardware design. 
Note that to get a patent it is necessary to present the so-called 'inven- 
tion formula', that is to indicate a prototype for the suggested device 
and the difference between these devices. 

Consider (as a typical example) the most popular and important 
algorithm of computing the scalar product of two vectors: 

(x,y) = xtyt + x 2 y2 + ■ ■ ■ + x n y n . (1) 

The universal version of (1) for any semiring A is obvious: 

(x, y) = (xi yi) © (x 2 y 2 ) © . . . © (x n y n )- ( 2 ) 

In the case A = M max this formula turns into the following one: 

(x,y) = max{xi + yi,x 2 + y 2 , . . . ,x n + y n }. (3) 

This calculation is standard for many optimization algorithms, so 
it is useful to construct a hardware unit for computing (3). There 
are many different devices (and patents) for computing (1) and every 
such device can be used as a prototype to construct a new device for 
computing (3) and even (2). Many processors for matrix multiplication 
and for other algorithms of linear algebra are based on computing scalar 
products and on the corresponding 'elementary' devices respectively, 
etc. 

There are some methods to make these new devices more universal 
than their prototypes. There is a modest collection of possible opera- 
tions for standard numerical semirings: max, min, and the usual arith- 
metic operations. So, it is easy to construct programmable hardware 
processors with variable basic operations. Using modern technologies 
it is possible to construct cheap special-purpose multiprocessor chips 
implementing reliable, thoroughly tested algorithms. The so-called sys- 
tolic processors are especially convenient for this purpose. A systolic 
array is a 'homogeneous' computing medium consisting of elementary 
processors, where the general scheme and processor connections are 
simple and regular. Every elementary processor pumps data in and out 
performing elementary operations in such a way that the corresponding 
data flow is kept up in the computing medium; there is an analogy with 



the blood circulation, hence the name 'systolic' (see, e.g., IS], fl32||). 

Of course, hardware implementations for important and popular 
basic algorithms can increase the speed of data processing. 
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9. Correspondence principle for software design 

Software implementations for universal semiring algorithms are not so 
efficient as hardware ones (with respect to the computation speed) but 
are much more flexible. Program modules can deal with abstract (and 
variable) operations and data types. Concrete values for these opera- 
tions and data types can be defined by the corresponding input data. 
In this case concrete operations and data types are generated by means 
of additional program modules. For programs written in this manner it 
is convenient to use special techniques of the so-called object oriented 



(and functional) design, see, e.g., pq] , [16|. Fortunately, powerful tools 
supporting the object-oriented software design have recently appeared 
including compilers for real and convenient programming languages 
(such as C ++ and Java). 

We propose a project to obtain an implementation of the corre- 
spondence principle approach to scientific calculations in the form of a 
powerful software system based on a collection of universal algorithms. 
This approach ensures working time reduction for programmers and 
users because of software unification. The arbitrary necessary accuracy 
and safety of numeric calculations can be ensured as well. 

The system contains several levels (including programmer and user 
levels) and many modules. 

Roughly speaking, it is divided into three parts. The first part con- 
tains modules that implement domain modules (finite representations 
of basic mathematical objects). The second part implements universal 
(invariant) calculation methods. The third part contains modules im- 
plementing model dependent algorithms. These modules may be used 
in user programs written in C ++ and Java. 

The following modules and algorithm implementations are in progress: 

— Domain modules: 

infinite precision integers; 
rational numbers; 
finite precision rational numbers; 
finite precision complex rational numbers; 
fixed- and floating-slash rational numbers; 
complex rational numbers; 

arbitrary precision floating-point real numbers; 
arbitrary precision complex numbers; 
p-adic numbers; 
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• interval numbers; 

• ring of polynomials over different rings; 

• idempotent semirings i?(max, min), i?(max, +), i?(min, +), 
interval idempotent semirings 

• and others. 
— Algorithms: 

• linear algebra; 

• numerical integration; 

• roots of polynomials; 

• spline interpolations and approximations; 

• rational and polynomial interpolations and approximations; 

• special functions calculation; 

• differential equations; 

• optimization and optimal control; 

• and others. 

This software system may be especially useful for designers of algo- 
rithms, software engineers, students and mathematicians. 
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